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ABSTRACT 

The  governing  equations  of  combustion  are  derived  from  the  standpoint 
of  continuum  mechanics;  and  an  indication  is  given  of  how  their  rational 
analysis  can  be  based  on  cctivation-energy  asymptotics. 
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Work  Unit  Number  3 (Applications  of  Mathematics) 


SIGNIFICArJCE:  AND  EXPLANATION 


A simi)le  but  careful  derivation  is  given  of  the  partial  differential 
equations  describing  the  flow  of  a chemically  reacting  mixture.  Simplify- 
ing assumptions  are  then  made  so  as  to  analyse  these  equations.  It  la 
finally  shown  that  the  largeness  of  the  so-called  activation  energy  in  many 
combustion  processes  enables  asymptotic  methods,  very  familiar  to  fluid 
dynamicists,  to  be  used  in  their  solution. 


The  r»s[x>n3 ibi 1 1 ty  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  arxl  not  with  the  authors  of  this  report. 


MATHEMATICAL  THEORY  OF  LAMINAR  COMBUSTION 


I.  GOVERNING  EQUATIONS 
>1.  D.  Bucluiiaster  and  G.  S.  S.  Ludford 

1 ■ Inf rodurt ion 

There  are  two  aspects  to  the  development  of  the  governing  equations  of  combustion: 
derivation  of  the  equations  of  a chemically  reacting  mixture  and  then  judicious  simplifica- 
tions to  render  these  equations  tractable  while  retaining  essential  characteristics.  A 
rigorous  derivation  requires  a long  apprenticeship  in  either  kinetic  theory  or  continuum 
mechanics.  (Indeed,  the  general  continuum  theory  of  reacting  mixtures  is  only  just  being 
perfected.)  We  choose  instead  the  path  of  a plausible  (but  potentially  rigorous)  deriva- 
tion, guided  by  experience  with  a single  fluid,  through  the  continuum  theory  of  a mixture 
of  fluids.  Such  a derivation  is  not  common  in  the  literature,  which  tends  to  run  to  ad  hoc 
argument  when  it  cannot  face  either  of  the  rigorous  ones.  In  particular,  the  inconsistent 
assumption  that  the  mixture  itself  is  a fluid  for  the  purpose  of  introducing  certain 
constitutive  relations  will  not  be  made. 

The  combustion  approximation,  characterizing  it  as  an  essentially  isobaric  process, 
is  a safe  simplification.  But  the  remaining  simpl if icatibns,  designed  as  they  are  solely 
to  make  the  equations  tractaule,  should  be  accepted  tentatively:  they  are  always  revoc- 
able in  the  light  of  faulty  predictions.  For  that  reason  they  should  be  explained  care- 
fully, but  are  usually  not.  Nevertheless,  whoever  is  primarily  interested  in  solving  non- 
trivial combustion  problems  (as  we  are)  can  have  the  same  confidence  in  the  final  equations 
as  is  normally  placed  in  the  Navier-Stokes  equations. 

These  final  equations  retain  most  of  the  complexity  of  a compressible,  heat-conducting, 
viscous  fluid  but  add  diffusion  of  the  species  and  source  terms  representing  the  chemical 
reaction.  This  complexity  has  usually  been  fought  with  irrational  approximation  and 
computers.  However,  combustion  processes,  by  their  very  nature,  tend  to  have  large 
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«ctivation  energies.  It  is  therefore  more  natural  for  a fluid  dynamicist  to  analyze 
the  problem  rationally  by  localizing  the  reaction,  for  example,  in  layers  (or  flames), 
which  Sec.  8 will  show  can  be  done  through  activation-energy  asymptotics.  Such  an  approach 
was  scarcely  mentioned  before  the  review  article  of  Williams  (1971),  which  contains  the 
first  suggestion  of  so  forming  the  basis  of  a mathematical  theory  of  combustion.  To  be 
sure,  the  classic  paper  of  Bush  & Fendell  (1970)  came  earlier,  but  Williams  mar)(ed  the 
start  of  a clear  path  through  the  complex  and  fascinating  problems  of  combustion,  at  least 
for  the  present  writers. 

Our  intention  then  is  to  arrive  as  quicisly  as  possible  at  a convincing  set  of  equa- 
tions and  a method  of  attaclc,  with  which  to  get  on  with  discussing  combustion  phenomena. 


p 

2.  Continuum  Approach  to  a Mixture  of  Reacting  Spec i^s 

The  mixture  has  density  u arxl  is  considered  to  be  made  up  of  N fluids  whose 
. separate  densities  are  CY^(i  • 1,2,...,N).  Here  the  Y.  are  mass  fractions  (or  concen- 

trations) , with 

(1)  I Y = 1 . 

i-1  * 

If  IS  the  velocity  of  the  ith  fluid  (or  species),  balance  of  mass  requires 

(2)  3(0Y.)/3t  + V • (dY.v.)  > 0.  , 

i - 1 -i  1 

where  is  the  rate  of  production  of  species  i (mass  per  unit  volume)  by  the  chemical 

reactions  so  that 

N 

(3)  I D.  =■  0 . 

i*l 

Summation  of  the  equations  over  all  species  therefore  yields  the  familiar  continuity 
equation 

(4)  3o/3t  + V • (Dv)  » 0 , 
where 

N 

' (5)  Y ■ I ’^iYi 

i-1  ^ 

is  the  mass-average  velocity  of  the  mixture.  Each  balance  equation  (2)  can  now  be  re-cast 
by  subtracting  Y^  times  the  overall  equation  (4)  and  introducing  the  diffusion  velocities 

(6)  V - V v ; 
this  gives 

(7)  0(3Y^/3t+y7Yj)  - p,  - 7 • (pY.Vj)  . 

The  moir.entum  balances  of  the  separate  species  will  not  be  written  down  since  they 
play  no  role  in  the  sequel.  But  Truesdell  (19fi5)  shows  that  if  they  are  added  the  result 
is  the  single-fluid  balance 

(8)  p ( 3y/3t*y -yv)  - y-r  + pij  , 

provided  no  mrxnentum  is  created  by  the  chemical  reactions.  Here  I is  the  sum  of  the 
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separate  stresses  plus  stresses  due  to  diffusion  of  the  species  and  g is  the  gravity 
force,  assumed  to  be  the  only  force  external  to  the  mixture.  Interactions  between  the 
species  [cf.  the  discussion  of  Pick's  law  (19)),  which  give  additional  external  forces 
on  each  individual  species,  sum  to  zero. 

We  shall  also  forego  writing  the  energy  balances  for  the  sei^rate  species,  but  Truesdell 
shows  their  sum  gives  the  single-fluid  balance 

(9)  o ( HU/it+vVU)  » I:Vv  + , 

provided  no  energy  is  created  by  the  chemical  reactions  and  the  work  of  the  interaction 

forces  is  negligible.  Here  U is  the  sum  of  the  separate  internal  energies  plus  kinetic 

energies  of  diffusion  while  q is  the  sum  of  the  separate  energy  fluxes  plus  fluxes  due 

to  diffusion.  Note  that  the  work  of  the  gravity  force  has  not  been  neglected:  it  is 
N 

I (pY.)gV  ”0  by  virtue  of  the  result 
i-1  ^ 

N 

(10)  I Y.V.  - 0 , 

i-1 

which  follows  from  the  definition  (6) . 

The  new  variables  o , V.,  E,  U and  q introduced  by  these  tialances  have  to  be 
1 **  1 « *• 

related  to  the  basic  variables  by  constitutive  equations.  Continuum  mechanics  does  not 
supply  Such  relations  but  rather  judges  them  for  consistency  with  certain  general  principles 
after  they  have  bean  proposed.  For  reacting  mixtures  such  judgements  are  still  being 
made,  but  the  linear  equations  we  shall  write  down  are  acceptable  (Bowen,  1976).  With  the 
exceptions  of  and  they  come  from  experience  with  a single  fluid.  The  coefficients 

(37)  which  will  appear  are,  for  the  moment,  to  be  considered  functions  of  the  basic 
variables. 

The  basic  variables  will  be  Y^,  v,  o and  p;  but  in  formulating  the  constitutive 
equations  it  is  convenient  to  introduce  the  temperature  T,  which  is  assumed  to  be  the 
same  for  all  species.  If  each  species  is  a perfect  gas  its  partial  pressure  is 

(11)  - R(pYj^)T/m^  , 
and  the  consequence 
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(12) 


p - R&T  I (Y,/m. ) 

i-1  i * 


of  Dalton's  law  determines  the  common  temperature.  Here  R is  the  gas  constant  and  m^ 
is  the  mass  of  a molecule  of  species  i.  The  separate  internal  energies  of  the  species 
are  then  functions  of  T alone,  which  may  be  written 
(13) 


U,  - h.  - p./(dY.) 
1111 


where  the  enthalpies  are 


(14) 


h.  - h + / c . (T)dT  . 

1 1 O 

T 

o o 

Here  h^  is  the  heat  of  formation  of  the  species  i at  some  standard  temperc.ture  T 

and  the  c^^  are  the  specific  heats  at  constant  pressure.  If  we  neglect  )cinetic  energies 

of  diffusion  for  being  nonlinear  in  the  V^,  the  internal  energy  of  the  mixture  is  now 


(15) 


0 - I Y.h  - p/o  . 

i-1  ^ 


Next  the  energy  flux  g is  specified  by  requiring  that  the  separate  energy  fluxes  be  due 
entirely  to  heat  conduction,  so  that 

N 


(16) 


q « IVT  - 0 y Y.h^V. 
‘ - **  1 1 *1 
i>l  ^ ^ 


wl.ere  1,  the  coefficients  of  thermal  conductivity  of  the  mixture,  is  the  sum  of  the 
species  coefficients.  Each  of  the  diffusion  fluxes  is  that  of  a single  fluid  moving 
relative  to  the  mixture  with  velocity  when  )(inetic  energy  and  deviation  of  the  stress 

tensor  (17)  from  neglected,  the  latter  anticipating  the  combustion  approximation 

(Sec.  5)  which  ma)(es  v^  small.  Such  a q neglects,  in  particular,  radiative  transfer 
and  the  Dufour  effect  (heat  flux  due  to  concentration  gradients) . 

The  separate  stresses  are  assumed  to  be  Newtonian,  i.e. 

(17)  - -(p.+2x^V-Vj/3) I ♦ ' 

where  I is  the  unit  tensor  and  bul)(  viscosity  has  been  neglected.  We  shall  ta)ce  the 
remaining  coefficients  amounts  to  supposing  that  the  intrinsic  viscosities 


of  the  species  are  all  equal.  (Otherwise  gradients  of  must  be  neglected  in  comparison 

with  those  of  v to  arrive  at  the  result  (18).)  Since  stresses  due  to  diffusion  and 

terms  V.Vy.  are  nonlinear,  these  sum  to  give 
■ 1-  1 


(18)  t - -(p+2icV.v/3)  r + «r(7v+(7v)  J 

by  virtue  of  the  result  (10) . 

The  thermomechanical  constitution  is  completed  by  equations  for  the  diffusion 
velocities  V^;  we  adopt  Onsaqer's  (1945)  generalization 


(.V.  ” - f u. .VY^ 
1-1  13-  j 


of  Fic){'s  (1855)  law,  malting  each  diffusive  mass  flux  a linear  combination  of  the  concentra- 
tion gradients.  Now  the  condition  (1)  and  the  result  (10)  imply  that 


N N N-1  r N ■) 


where  the  remaining  greulients  are  independent  of  each  other.  It  follows  that  the 

coefficients  u, . must  satisfy 
il 


I U - = I I'm  for  3 “ 1,2,...,  N-1  . 


i«l  i=l 


An  equivalent  of  the  law  (19),  in  which  the  gradients  are  expressed  as  linear  combina- 
tions of  the  diffusion  velocities,  can  be  derived  from  the  separate  momentum  balances  by 
malting  assumptions  about  the  interaction  forces  between  the  species  (Williams  1965,  p.  416)  . 
An  implicit  new  assumption  in  such  arguments  is  tliat  the  accelerations  3v^/3t  are 
negligible.  They  only  occur  in  unsteady  problems  and  then  produce  a finite  velocity  of 
propagation  of  diffusion  effects  (Muller  1977);  as  such,  they  appear  to  be  of  no  great 
importance  in  combustion.  The  constitutive  equation  (19)  is  the  only  non-chemical  one 
which  does  not  come  from  experience  with  a single  fluid.  Note  that  the  rate  of  working 
of  these  interaction  forces  is  nonlinear  in  so  that  we  are  justified  in  using 

equation  (9) . 

Finally  we  come  to  the  o^,  which  require  a discussion  of  chemical  reactions. 
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3.  The  Arrhenius  F'actor  (1889) 


For  slm^licity  we  shall  first  consider  one-step  combustion,  where  a sinqle  unopposed 
chemica,^  reaction  is  involved.  If  is  the  number  of  molecules  of  species  i per  unit 

volume,  then  (by  definition) 

*21 ) r . " m . N . 

1 1 i 

The  chemical  reaction  may  be  described  in  terms  of  the  m.  as  the  mass  balance 

1 


(22) 


IN  M 

y v.m.  - y Am. 
i-1  ^ ' i-1  ^ * 


where  coefficients  v^,  A^  are  the  stoichiometric  integers,  each  a non-negative  integer, 
v^  is  zero  if  the  species  i is  not  a reactant,  while  A^  is  zero  if  it  is  not  a product. 
Both  are  zero  for  an  inert.  Since  these  coefficients  represent  the  relative  numbers  of 
molecules  consumed  or  produced  by  the  reaction  we  may  write 


(23) 


N.  = (X  -u,)i, 
1 i 1 


where  u measures  the  rate  at  which  the  reaction  is  proceeding.  Clearly  the  requirement 
(3>  is  now  a consequence  of  the  definition  (21)  and  the  mass  balance  (22). 

For  example,  consider  the  exothermic  reaction 


(24) 


2N0  * Cl  ■*  2N0Cf 
2 


by  which  nitrosyl  chloride  is  formed  from  nitric  oxide  and  chlorine.  If  the  three  species 
are  numbered 


1:N0  2!CI^  3:N0Cl 


then  the  stoichiometric  integers  are 


''l  " '*2  ■ '*3  “ - 0,  >2  > 0,  Xj  - 2 , 


and 


” -2u,  Nj  “ -u),  Nj  • 2u 

are  the  rates  of  consumption  (negative)  and  production  (positive)  of  molecules.  Note  that, 
when  single  molecules  are  used,  the  reaction  (24)  is  written 

NO  ♦ NO  + Cfj  Noce  + NOC£  . 
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It  is  therefore  u)  which  must  be  relatfxi  to  the  basic  variables;  a common  assumption 
is  that  it  is  proportional  to  the  number  density  (concentration)  oYVro^  of  each  reactant, 
where  the  species  j is  counted  times  as  a reactant  (cf.  the  single-molecule  form 

of  the  reaction  (24) ) . Thus 


(25) 


<D 


N 

w n 


j-i 


N ( Y . 1 ^ N 

ItO  1 r where  u “ 7 u. 

j-irjl  A ^ 


and  the  factor  k is  supposed  to  depend  only  on  the  temperature.  (Since  experimental 
determinations  of  k are  done  under  isothermal  conditions,  it  is  called  the  rate  constant; 
but  we  shall  avoid  the  term  since  the  temperature  dependence  of  k is  an  essential  feature 
in  our  analysis.)  The  Arrhenius  law 


(26) 


, -E/RT 

k = BT  e 


is  invariably  taken  for  its  variation  with  temperature,  where  B,  a (*  1/2  theoretically), 

E are  constants.  The  latter  is  known  as  the  activation  energy  since  E/R  is  roughly 

the  temperature  below  which  k is  small.  In  short,  we  now  have 

(27)  D,  ” u.u 

i 1 

where 


(28) 


U . 
1 


,,  , , a -E/RT  V V I 1 

mj(X^-vJ  and  u » BT  e p ^ 1 — ^ 


1 i i 


j»l  \"j 


More  complex  reactions  are  first  split  up  into  elementary  reaction  steps,  which  may 
contain  reversals  of  each  other  (opposing  reactions) . For  each  elementary  reaction  a 
formula  (27)  holds  and  for  the  complete  reaction  is  the  sum  of  all  such  terms.  (A 

pair  of  opposing  reactions  may  be  in  equilibrium,  i.e.  the  corresponding  terms  cancel  and 
thereby  furnish  an  algebraic  relation  between  mass  fractions.)  The  central  question  of 
reaction  kinetics  is  to  determine  these  elementary  steps  and  their  parameters  B,  a,  E. 

In  practice  a combustion  reaction  may  be  a complex  network  of  elementary  steps. 
Attempts  are  therefore  made  to  model  the  overall  reaction  with  a smaller  number  of  terms 
of  the  form  (27),  preferably  one.  The  exponents  v^  are  then  released  from  being  the 
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stoichiometric  coefficients  of  the  overall  reaction  and  v is  not  necessarily  their  sum. 

On  the  other  hand,  the  u . retain  their  meaninq  since  they  still  represent  the  proportions 


P 

• in  which  the  specli^s  take  part  in  the  reaction.  The  free  constants  B,  a,  E,  v,  are 

determined  experimentally;  of  course  the  need  no  longer  be  integers. 


I 

I 


4.  Differential  Dif fusion.  Equality  of  Specific  Heats  and  of  Molecular  Masses 

Once  the  various  parameters  (constants  and  functions)  which  appear  in  the  constitutive 

equations  (12),  (15),  (16),  (18),  (19),  (27)  are  given,  there  are  N 6 equations  (1), 

(4),  (7),  (8),  (9)  for  the  N + 5 unlcnowns  Y.,  v,  0,  p.  Equation  (1)  is  consistent 

N ‘ 

with  the  remainder  since  they  conserve  ) Y.  whenever  the  relations  (20)  hold  between 

i”l  ^ 

the  I'jj'  initial  or  boundary  conditions  will  then  give  it  the  value  1.  The  system  is 
clearly  very  complicated  as  it  stands,  and  must  be  simplified  if  it  is  to  yield  to  detailed 
analysis. 

The  law  (19)  leads  to  terras 

N 

y v.(ii,.vy.) 
j=l  ' ’ 


(29) 


in  the  mass  balances  (7) . To  simplify  these  terms  while  retaining  different  diffusion 

properties  for  the  species,  it  is  tempting  to  set 

(30)  u . . 0 for  i ^ j 

13 

and  allow  the  to  change  with  i (ESnmons  1971)  . However  the  conditions  (20)  would 


be  violated:  they  require 

(31) 


Wii  ’ ” l](say)  for  all  i . 


The  conclusions  (30)  and  (31)  are  reached  quite  generally  by  Williams'  theory  when  the  so- 
called  binary  diffusion  coefficients  are  equal,  i.e.  whenever  the  interactions  between 
pairs  of  species  are  identical  so  far  as  momentum  is  concerned. 

However,  effectively  the  same  goal  can  be  reached  by  taicing  the  values  (30)  for  all  i 
except  the  last,  and 

(32)  for  j # N 

with  the  different.  The  conditions  (20)  are  now  satisfied  and  the  terms  (29) 

simplify  to 

(33) 


V-  I * H . 


Only  the  Nth  equation  remains  complicated  knit,  since  Y can  be  otherwise  calculated 
N-l 

from  1 - y Y , that  does  not  matter.  Such  a scheme  can  be  justified  from  Williams' 
i-1  ‘ 
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formulas  whon  the  first  N-X  species  are  scarce  comj>arod  to  the  last,  an  extromf*  often 


met  in  coirtustion  where  , for  example,  reactants  are  diluted  in  an  inert  species  so  as  to 
prevent  explosions. 

In  short,  we  shall  generally  taxe  the  values  (30)  and  (31)  but,  when  differential 

diffusion  cannot  be  ignored,  the  u , (i  = 1,2,...,N-1)  will  be  given  different  values. 

11 

The  results  will  then  apply  to  a mixture  of  reactants  and  products  highly  diluted  in  an 

inert.  Note  that  a common  value  of  the  u. , can  still  depend  on  the  Y.  and,  in 

ii  1 

particular,  be  a different  constant  when  different  Y.  vanish.  Such  is  the  case  in  the 

1 

burning  of  a fuel  droplet  treated  by  Kassoy  & Williams  (1968).  They  use  different  diffu- 
sion coefficients  for  the  fuel  and  oxidant  equations  but  do  not  apply  the  equations  in 
the  same  place:  the  combustion  field  is  divided  into  two  parts,  in  each  of  which  there 
is  either  no  fuel  or  no  oxidizer.  There  is  no  question  of  an  abundant  inert. 

We  come  now  to  the  energy  balance  (9) . First  note  that  the  individual  disappear 

from  U in  favor  of 


(34) 

Thus  we  have  immediately 

(35) 


c (T, y.) 
P 1 


N 


2 

i»l 


Y 


. c , . 

1 pi 


N 

U*  y Y.h?+I-  p/o 
I 1 


where 


(36) 


NT  T 

I(T,Y,)  - y Y,  f c dT  * I c dT 
1 . 1 •'  pi  P 

T T 


the  last  integral  being  taken  with  the  Y^  fixed,  i.e.  for  constant  composition.  Adoption 
of  ensures  that  the  same  happens  to  q:  we  have 


(37) 


N N N ; T \ 

y Y.h.V.  - 0 y Y.h°V.  - M y / c .dT  VY. 

i i-i  i i-i  p'  pi  ' I 

1-1  i-l  1-1  1 _o  / 

T 

N N I 

- 0 y Y.h°V.  - u 71-  y Y.c  .7T 

i-l  ‘ ' 1-1  ‘ I 


0 y Y,h°v.  ♦ y(c  7T-VI)  . 
iil  i P-  - 
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Next  note  the  form  which  the  energy  balance  now  takes.  The  total  contribution  to 


the  equation  from  the  heats  of  formation  h^,  when  placed  on  the  right-hand  side,  is 

N N 

(38)  - I h°|p(3Y./at-v-VY.)+V.(DY  v.)l  --I  h°o . 

i"!  ^ * i«l  ‘ 

accordinq  to  t)>e  mass  balances  (7)  . Also  the  total  contribution  from  heat  conduction  and 

terms  contalninq  c^,  when  placed  on  the  left-hand  side,  is 

(19)  p(3I/at+v.7I)  ♦ v-|u(c  7T-Vl)-iVTI  . 

p,  . 

The  remaining  terras,  when  placed  on  the  right-hand  side,  contribute 

(40)  3p/3t  + v-7p  - 2ic(7-v)^/3  * < (7v+ (7v)^l  : 7v 

according  to  the  continuity  equation  (4) . In  the  combustion  approximation  discussed  latei 
the  velocity  is  small.  Implying  also  that  the  pressure  is  effectively  a spatial  constant, 
so  tliat  these  last  terms  are  ignored  except  for  3p/9t,  which  is  a function  of  t alone. 
It  is  tempting  to  ta)(e  the  Lewis  number 

(41)  i » 1/uc 

P 

equal  to  Ij  then  the  terms  (39)  contain  I only.  However,  if  the  c , are  not  all  equal, 

pi 

c is  a linear  function  (34)  of  tiie  Y,  which  u/1  can  only  equal  for  very  special 
P i 

mixtures  (if  at  all).  Only  when  the  c^^  are  all  equal  does  c^  lose  its  dependence  on 
the  composition  of  the  mixture,  and  then  it  is  unnecessary  to  ma)ce  any  assumption  about 
the  Lewis  number  as  we  shall  see  next. 

We  are  therefore  faced  with  tailing  all  the  c , equal,  when  c becomes  a function 

pi  p 

of  T alone: 

(42)  c (T)  “ c , for  all  1 . 

P pi 


Since  7l  is  than  c^7t  the  energy  equation  reads 


>(3l/at+y-7I)  - 7,l(l/Cp)7Il  * I *’i‘‘i]“  * 


I(T)  - / c dT 

;«> 


is  a function  of  T only  and  we  have  used  (27).  So  far  as  specific  enthalpy  is  concerned 
the  species  are  distinguished  only  by  their  heats  of  formation;  thus  equation  (43)  is 
that  of  a single  fluid  in  which  differences  in  these  heats  are  liberated  as  reactants 


change  into  products.  [It  sliould  be  noted  that  the  same  equation  results  under  the  assump- 
tion of  equal  specific  heats  even  when  there  is  differential  diffusion.)  The  species 
equations  (7)  talie  a similar  form 


(45) 


D(3V./3t+v-VY,)  - V-{uVY.)  - u.u  . 
1 ' • 1 - ' 1 1 


Another  complication  is  the  presence  of  the  mass  fractions  in  the  aquation  of  state  (12). 
They  effectively  disappear  if  the  first  N-1  species  are  scarce  compared  to  the  last; 
the  same  simplification  can  only  be  achieved  in  general  when  all  the  molecular  masses  are 
the  same,  i.e. 


(46) 


m^  » m (say)  for  all  i . 


Then  the  law  (12)  becomes 

(47)  p . RpT/m 

so  that  the  mixture  believes  lilte  a perfect  gas.  Such  is  the  case  when  the  reaction  changes 
a single  molecule  into  a single  molecule,  but  otherwise  the  condition  is  never  met  exactly 
in  practice. 

We  may  expect  our  equations  to  govern  the  behavior  of  a mixture  for  which  c , and 

pi 

m do  not  vary  much  from  species  to  species.  In  practic  c and  m would  be  given 

N r N l-i 

intermediate  values,  e.g.  I Y c and  I Y./m  I for  a characteristic  set  of  Y. . 

i-1  P l-i-l  ^ ^ 

More  serious  is  the  consequence  (31)  of  the  assumption  (30).  Differential  diffusion  is 


linown  to  be  responsible  for  important  effects  in  biological  systems,  in  particular  pattern 
forsiation  (uthmer  1977);  and  is  believed  to  produce  instability  leading  to  the  cellular 
structure  of  flames  (Marlistein  1964,  p.  75).  Luc)illy  the  latter  question  can  be  tacliled 
with  the  alternatives  (32)  and  (33).  Either  there  is  an  abundant  inert  or  all  but  the  single 
product  species  (substituting  for  the  inert)  are  scarce  near  the  flame,  whose  vicinity 
alone  need  be  considered. 
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The  coefficients 


(461  «■,  i,  u.  c 

P 

are  functions  of  the  basic  variables  Y^,  0 and  p.  Further  sinplifications  will  result 
froai  takin9  then  constant/  but  we  shall  delay  writing  down  the  equations  until  those  are 
stade  diSMnsionless  later. 
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5.  The  Combustion  Approximation  and  Its  Constant-Dens ity  Approximation 

While  the  simplifications  introduced  so  far  have  the  nature  of  idealizations  which 
are  approached  more  or  less  by  actual  mixtures,  the  basic  approximation  of  the  theory 
characterizes  combustion  as  a low-speed  phenomenon  in  a chemically  reacting  compressible 
mixture  for  which  gravity  is  unimportant.  Mote  precisely,  the  kinetic  energy  is  small 
compared  to  the  thermal  energy  so  that  a representative  Mach  number  is  small. 

An  immediate  consequence  is  that  spatial  variations  in  pressure  are  also  small,  so 
that  we  may  set 

(49)  p - p^(t) 

everywhere  except  in  the  momentum  equation  (8) . (The  notation  indicates  that  the  pressure 
can  be  controlled,  for  example  in  the  surroundings.)  This  conclusion  is  reached  from  the 
momentum  equation,  which  shows  that  fractional  changes  in  pressure  at  any  time  are  propor- 
tional to  the  Mach  number  (both  'hrough  3v/3t  and  the  viscous  terms).  The  pressure  can 
still  vary  in  time,  but  a subsidiary  feature  of  the  combustion  approximation  is  that, 

when  it  does,  the  variations  are  small  ones  about  some  constant  value  p . 

c 

Buoyancy  v#ould  appear  to  be  an  important  factor  in  certain  combustion  problems.  If 
so,  pg  should  be  retained  in  the  momentum  equation  and  the  conclusion  (49)  cannot  be 
drawn.  Such  a problem  is  the  burning  of  a fuel  drop  but,  nevertheless,  remarkably  accurate 
predictions  of  combustion  characteristics,  as  opiosed  to  the  flow  field  (which  is  not  of 
primary  interest),  result  from  neglect  of  gravity.  (Cf.  Williams  (1965,  p.  48).)  The  reason 
may  lie  in  the  absence  of  gravity  from  the  energy  equation  (9)  and  in  an  overall 
insensitivity  of  its  solutions  to  inaccurate  determinations  of  the  convection  ov. 

It  is  tempting  to  think  that  the  momentum  equation  can  now  be  discarded  because  its 
only  role  is  the  calculation  of  these  small  "f lovf-induced"  variations  in  pressure.  But 
that  would  leave  N ♦ 2 equations  for  the  N ♦ 4 unknowns  Y^,  y,  o.  If  these  variations 
must  be  eliminated,  the  momentum  equation  should  be  replaced  by  Its  own  curl  (Helmholtz 
equation),  which  provides  two  Independent  scalar  equations.  Wo  prefer  to  leave  it  intact 
and  interpret  p as  the  small,  flow-induced  variation  in  pressure  from  its  value  p^. 


-15- 


r 

If  the  reaction  rate  u is  zero  (as  it  is  in  whole  regions  for  the  asymptotic 
approach  we  shall  describe  in  Sec.  8),  the  five  equations  (4),  (8),  (43)  for  the  five 
unknowns  v,  p,  o become  uncoupled  from  the  species  equations.  They  can  be  recognized  as 
the  equations  of  motion  of  a viscous,  heat-conducting,  compressible  fluid,  for  which  few 
exact  solutions  are  known.  To  be  sure,  the  combustion  approximation  has  been  applied  but 

such  modification  leaves  them  only  slightly  more  tractable.  A non-zero  u can  only  make  ^ 

analytical  solution  more  difficult,  which  accounts  for  the  heavy  reliance  on  numerical 
methods  in  the  combustion  literature. 

Most  progress  can  be  made  for  steady  conditions,  since  only  the  mass  flux  pv  is 
needed  to  make  the  energy  equation  (43)  and  species  equations  (45)  a determinate  set  for 
T and  Y^.  For  certain  geometries  the  mass  flux  is  obvious  a priori  while  for  others 
it  can  be  approximated  locally.  The  most  important  example  of  the  former  is  a uniform 
flux  in  the  x-direction,  for  which 

(50)  pv  • Mi,  M const. 

It  can  be  used  when  T varies  with  x alone,  but  otherwise  needs  justification.  If 
T also  varies  radially  from  the  x-axis,  then  so  do  p • p^m/RT  and  the  velocity.  While 

there  are  such  unidirectional  velocity  fields  satisfying  the  momentum  equation  these  are  • 

not  amongst  them  (as  is  easily  shown)  . Burke  t Schumann  '(1928),  in  tackling  what  may  be 

i 

regarded  as  the  first  combustion  problem,  ignored  this  inconsistency . 

Justification  comes  when  the  variations  in  temperature  are  due  entirely  to  the  combus-  | 

tion  and  the  heat  release  is  small  compared  to  the  thermal  energy  of  the  mixture.  Since  ’ 

the  density  variations  are  correspondin<jly  small,  the  velocity  field  is  that  of  a constant- 
density  fluid  to  leading  order.  The  undisturbed  values  of  P and  pv  may  be  used  in 
i equations  (43)  and  (45)  to  determine  the  leading  terms  in  the  mass  fractions  and  the 

^ perturbation  of  the  temperature.  In  short,  Burke  t Schumann  introduced  the  constant- 

^ density  approximation,  a rational  scheme  valid  even  for  unsteady  problems.  No  attempt 

has  been  made  to  take  it  any  further  than  they  did. 

. \ 

4 


I 
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6.  Constant  Properties.  Dimensionless  Variablea. 

The  coefficients  (48)  will  now  be  ta)<en  constant. 

The  startinq  point  for  makinq  the  variables  dimensionless  is  the  unit  of  tunperature, 
which  comes  from  the  factor 

N N 


(51) 


■ I h°li  -ml  (say) 

i-1  i»l 


on  the  right-hand  side  of  the  energy  equation  (43) . Q is  the  excess  of  the  total  heat  of 

formation  of  the  reactants  over  that  of  the  products  per  unit  mass  of  the  reactants  (or 

products),  i.e.  the  heat  per  unit  mass  released  by  the  reaction  at  the  standard  temperature 

T°.  [When  all  the  c . (T)  are  equal,  as  hare,  it  is  also  the  heat  released  at  any 
pi 

temperature.)  Combustion  is  only  concerned  with  exothermic  reactions  such  as  (24)  and 
then  Q(>  0)  is  called  the  (specific)  heat  of  combustion.  The  unit  of  temperature  is  ta)(en 
to  be  0/Cp<  so  that  the  perfect  gas  law  becomes 
(52)  T " 

if  pressure  is  referred  to  p^  and  density  to  " CpP^m/QR. 

The  pressure  p will  be  varied  from  0 to  «,  so  that  the  choice  of  density  unit 
c 

ensures  that  o stays  finite.  The  same  should  hold  for  the  mass  flux  but  there  is  no  such 

combination  of  existing  par^uneters  (excluding  the  inappropriate  rate  constant  B)  . We 

therefore  take  a representative  mass  flux  M and  refer  the  velocity  to  N/o  . A diffusion 

c 

2 

length  X/c  M and  a diffusion  time  o X/c  M can  then  be  defined  as  units,  in  recogni- 
P c P 

tion  that  these  are  the  scales  for  spreading  the  heat  generated  by  the  combustion.  Finally 

2 

pressure  variations  are  referred  to  H /p  . 

c 

In  these  units  the  governing  equations  become 


(53) 

(54) 

(55) 

(56) 


3o/3t  ♦ V • (pv)  • 0 , 


p(3Yj/3t+v.7Y^)  - L'^V^Yj  - a^n  , 


0 ( 3y/9ttv 'Vv)  “ -7p  + P(V  v+V(7.v)/3)  , 


DOT/3t+V'7T)  - 7 T - 0 ♦ B3p  /3t 
• c 
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where  the  dimensionless  rate  of  heat  release  per  unit  volume  is 

(57)  a - rr  Y A = Dm'^  . 

j-1  ^ 

Here  P » c k/A  is  the  Prandtl  number,  L has  the  definition  (41),  the  quantities  charac- 
P 

terizinq  the  chemical  reaction  are 

(58)  a.  - u,/vm  - (>,-v,)/v,  0 - Ec  /Vf>  , 

1 i i 1 p 

1*V  Ol“V  “\>“V  01  V 

(59)  D(T)  - (Bvm  AQ  c R p )T  p , 

P c 

and  B • (yDA  where  y « l/(l-mc  /B)  is  the  ratio  of  specific  heats. 

P 

When  the  geometry  provides  a characteristic  length  a the  nondimensionalization  is 

2 2 2 2 

usually  done  differently.  A/p  c a,  p c a /A  and  A /p  c a are  ta)cen  as  units  of 

cpcp  cp 

velocity,  time  and  pressure  variation,  respectively,  to  give  the  same  equtions  (53-56) 
except  that  now  M is  missing  from  the  definition  (57)  of  (2  and 

1-v  -l.u-v  v-a+l„-v  2-v,_a  v 

(60)  D(T)  - (Bum  Afi  c Rap)Tp. 

p c 

In  effect  M has  been  replaced  by  the  new  unit  of  mass  flux  A/c^a. 

For  any  species  which  is  solely  a reactant/product  is  negative/positive;  also 


I “i  - 0 , 


since  the  sum  to  zero  according  to  the  mass  balance  (22).  The  variations  in  D(T) 

are  unimportant  and  its  value  at  some  intermediate  temperature  may  be  ta)ien;  in  our  analysis 
only  its  value  at  the  flame  sheet  will  appear.  D may  then  be  called  the  Dam)cohler  number, 
representing  the  pressure  p^.  Finally,  0 is  the  dimensionless  activation  energy. 

When  the  reactants  and  products  are  distinct,  so  that  we  may  talte 
• 0,  v^  # 0 for  i - 1,2, ...,n  (say)  , 


Vj  ■ 0 


for  1 " n+l,n+2, . . . ,N  , 


the  equations  tor  the  products  and  inerts  can  be  solved  after  the  rest.  Thus  the  first  n 
of  the  species  equations  (54)  contain  only  the  reactant  mass  fractions,  as  does  the  energy 
aquation  (56);  together  with  the  equations  of  continuity  (53)  and  sximentum  (55)  they  form 
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n ♦ 5 equations  for  the  n + 5 unknowns  Y,,Y,,...,Y  , v,  p and  p.  The  coefficients  on 

1 2 n ■ 


the  right-hand  sides  of  these  reactant  equations  have  the  property 


y a.»-  y v./u  ■-! 
..  i 1 

i«l  I'l 


Finally  we  mention  the  modifications  when  differential  diffusion  is  considered:  the 

N-1 

inert  species  equation  i • N is  dropped  in  favor  of  Y = I - y Y.  and  L is  allowed 

i»l 

to  vary  among  the  remaining  species. 


7 . Shvab-Zeldovich  Variables 


When  the  Lewis  number  (41)  is  1,  the  convection-diffusion  operators  on  the  left-hand 
sides  of  equations  (54,56)  are  identical,  so  that 

(64)  (o3/3t+Ov.V-V^) (V^-a.T)  -0  for  all  i . 

The  reaction  terms  have  been  eliminated  to  show  that  the  Shvab-Zeldovich  variables 

(65)  ’'^i  * ’^i  ■ “i'’’ 

and  purely  converted  and  diffused  as  in  an  inert  mixture.  The  can  often  be  determined 

first,  leaving  only  the  temperature  equation  (56) , with  the  Y^  replaced  by  Y^  + q^t  in 
n,  to  be  solved.  There  is  therefore  considerable  advantage  in  setting 

(66)  Z,  = 1 . 

However,  that  is  not  the  reason  for  ta)ting  the  value  (66).  The  asymptotics  of  large 
activation  energy  require  SI  to  be  zero  everywhere  except  in  certain  flame  sheets,  so 
that  there  are  only  convection-diffusion  equations  to  solve  there  anyway.  Inside  the 
flames  the  convection  terms  are  negligible  and  the  operators  are  once  more  identical  (except 
for  the  factor  L ^):  local  Shvab-Zeldovich  variables  can  be  formed  for  arbitrary  Lewis 
number.  The  excuse  for  the  value  (66)  is  solely  to  eliminate  an  aw)cward  symbol  and,  when- 
ever significantly  different  phenomena  can  be  obtained  (as  in  unsteady  problems),  a general 
value  should  be  adopted. 

These  convection-diffusion  equations  have  a particularly  simple  form  when  conditions 
are  steady.  If  the  (dimensionless)  energy  flux 

(67)  E • PTV  - Vt 
is  introduced,  then 

(68)  7-E  - 0 

in  view  of  the  continuity  equation  (53).  Similarly  the  (dimensionless)  species  fluxes 

(69)  J,  - py.v^  - 0Y,(v+VJ  - pvY,  - i'^VY^ 

•i  i'i  i - -i  -1  - i 

satisfy 

(70)  7-Jj  - 0 . 
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Equation  (68)  even  holds  when  conditions  are  unsteady,  provided  is  constant,  by 

virtue  of  equation  (52).  (Of  course  these  are  nothing  more  than  conservation  equations.) 

In  one  dimension,  integration  is  immediate!  and,  even  when  reaction  is  present,  the  species 
fluxes  are  useful  as  intermediate  variables. 


8.  Activation-Energy  Asymptotics 


The  reaction  (24)  has  an  activation  temperature  E/R  of  about  2,500‘’K  for  the  temperature 
range  273-523°K  (Kontrat'ev  1972)  so  that  exp(-9/T)  is  at  most  8 » 10  To  be  sure, 
the  reaction  is  exceptional  in  that  stoichiometric  integers  are  nearly  always  0 or  1 

and  sum  to  at  most  2 (corresponding  to  unimolecular  or  bimolecular  rather  than  termolecular 
reactions);  it  was  chosen  to  exhibit  the  square  of  a concentration  in  w.  But  it  is  not 
exceptional  in  having  a relatively  large  activation  temperature;  the  limit  0 -►  »>  is 
particularly  appropriate  for  many  combustion  reactions. 

To  see  an  immediate  effect  of  letting  0 -►  »>  it  is  convenient  to  write  the  reaction 
term  (57)  in  the  form 

v^  v V 

(71)  a = A'y/y^^ V.,^exp  e(l/T.-l/T),  A’  = A exp(-e/T.) 

12  N * 

where  T„  is  a characteristic  temperature,  as  yet  unspecified.  [It  will  be  the  temperature 
of  the  flame  sheet  to  which  the  chemical  reaction  is  confined  in  the  limit.]  We  now  assume 
A'  is  at  most  algebraically  large  in  0,  i.e.  that  all  the  exponential  growth  of  A has 
been  removed  by  the  factor  exp(-0/T^) . The  assumption  is  justified  if,  for  example,  it 
leads  to  a consistent  determination  of  D as  a function  of  the  burning  rate  M,  which 
is  fixed  as  0 increases. 

TVo  possibilities  now  arise.  The  frozen  limit  is  obtained  wherever  all  reactants  are 
present  but  T is  smaller  than  T^^ . The  term  n is  exponentially  small  so  that  to  all 
orders  in  1/0  there  is  negligible  reaction.  The  equilibrium  limit  obtains  v;herever  T 
is  larger  than  T,.  Thus 

(72)  - 0 for  some  i , 

to  all  orders  in  1/0,  since  otherwise  S)  would  )3e  exponentially  large.  Again  there  is 
negligible  reaction,  but  now  because  one  of  the  reactants  has  been  usecTup.  Of  course  the 
remaining  reactants  are  still  available  and  they  behave  as  if  the  chemistry  were  frozen. 

Only  where  T is  approximately  equal  to  T^,  more  precisely  where  T - T^  ” 0(0  ), 

does  a contribute;  so  that  the  reaction  is  confined  to  thin  regions,  iinown  as  flame  sheets. 


and  is  now  identified  as  the  flame  temperature.  On  either  side  the  reaction  iruay  be 

frozen  or  there  may  be  equilibrium.  If  the  reaction  is  frozen  on  both  sides  of  the  flame 
sheet,  the  ccxnbustion  is  not  complete. 

It  is  the  flame  sheet  which  determines  AMT^^)  (and  hence  as  a function  of  M)  ; 

for  each  M there  is  just  one  giving  a structure  which  can  be  matched  to  the  limit 

states  on  either  side.  Note  how  the  temperature  dependence  of  D has  dropped  out:  it  is 
treated  as  a constant  in  the  analysis.  The  general  procedure  in  any  particular  problem 
is  now  clear.  The  c«nbustion  field  is  divided  into  two  parts  by  a flame  sheet  and  either 
equilibrium  or  frozen  reaction  is  assigned  in  each.  (The  appropriate  assignment  is  rarely 
obvious  in  advance,  so  that  the  consequences  of  all  possible  combinations  must  usually  be 
worked  out.J  The  two  states  are  then  determined  as  far  as  possible,  their  final  determina- 
tion being  a result  of  matching  with  the  structure  of  the  flame  sheet.  Different  structures 
are  obtained  for  equilibrium  on  both  sides,  one  side  or  neither.  In  certain  circumstances 
the  flame  may  be  at  a boundary  (including  infinity) . Then  there  is  only  one  state  to  be 
determined  and  new  structures  arise. 

There  are  other  possibilities  for  the  asymptotics.  M and  D may  both  be  given  and 
some  other  parameter  is  to  be  determined;  then  D’  may  be  exjx^nentially  large.  Slow 
variations  of  otherwise  steady  states  (determined  as  above)  can  occur  on  a time  scale  0(0). 
In  so-called  thermal  runaway  the  time  scale  may  even  be  exponential  in  0. 


r 
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